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I. INTRODUCTION 



Ferroelectrics and the closely related high dielectric con- 
stant materials are important in many areas of modern tech- 
nology including memory, sensor and electronic applications, 
and are of fundamental scientific interest!! The ferroelec- 
tric phase change belongs to a class of structural transitions, 
generally termed ferrodistortive, triggered by zone-center soft 
modes of lattice motion. Characteristically the ferroelectric 
transition involves the condensation of an optically active lat- 
tice mode which causes the appearance of long-range polar 
order and the breaking of the inversion symmetry of the 'high- 
temperature' prototype lattice. One important issue is the fer- 
roelectric quantum critical point, i.e. the physics occurring 
when, by varying a control parameter r (applied pressure or 
change of chemical composition), the transition temperature 
of a ferroelectric is driven to zero. Although a quantum criti- 
cal point occurs at T = 0, the fluctuations associated with the 
critical point may control behavior over a range of tempera- 
ture and pressure. For example K. A. Mullera coined the term 
"quantum paraelectric" to describe materials in which ferro- 
electric ordering is prevented by quantal fluctuations. A mate- 
rial just on the disordered side of the T = ferro-paraelectric 
transition is therefore an example of a "quantum paraelectric". 
The important feature of quantal (T = 0) phase transitions is 
that temporal fluctuations must be treated on the same footing 
as thermal onesBo. This raises the effective dimensionality 
and makes the critical behavior more mean-field-like but with 
temperature dependence controlled by 'dangerous irrelevant 
operators '0. 

Ferroelectric transitions may be described by bosonic field 
theories with undamped dynamics (if the effect of free car- 
riers may be neglected) and complicated dispersions arising 
from the long range of the dipolar interaction. Quantum criti- 
cal phenomena associated with undamped bosonic field the- 
ories with short ranged interactions have been extensively 
studiedo. The-effect of long range (dipole) forces was-,studied 
by Rechesteru and by Khmel'nitskii and ShneersonEI within 
an approximation equivalent to the self-consistent one loop 
approximation of MoriyaH Aharony and Fishem studied the 



classical ferroelectric transition and found that anisotropics 
associated with the dipolar interaction led to a new univer- 
sality class. In this paper we reexamine the issue in light of 
recent developments in the theory of quantum critical phe- 
nomena. We formulate a realistic action for the ferroelectric 
soft modes, show how estimates of the parameters may be ob- 
tained from ab initio calculations and study quantitatively the 
consequences of the dipolar-induced anisotropies. Our results 
agree in essentials with those of Rechester and Khmel'nitskii 
and Shneerson, but we obtain a more detailed and quantitative 
picture of the phase boundary, of the effect of anisotropy, and 
of the logarithmic corrections arising at the marginal dimen- 
sionality, which lead to an evolution of the anisotropy as the 
ordered phase is approached. 



II. ORDER PARAMETER AND ACTION 

The order parameter of our theory is the local polarization 
4>. Taking into account the effective dipole charges e* of the 
soft modes, <fi can be formally written 
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Here the index i runs through the atoms of the unit cell of the 
prototype perovskite lattice with stoic hiometry ABO3; are 
the vector displacements of each atom. 

We now write a Ginsburg - Landau action describing quan- 
tal and thermal fluctuations of <p(x,t). The crucial point is 
that because <p corresponds to a dipole fluctuation it generates 
electric fields which lead to a long-range interaction. We have 
(in space and imaginary time) 
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Here a is the lattice constant, c is the speed of the phonons 
in the softest direction, and Eq = hc(ir/a) is the typical 
energy scale of ferroelectric fluctuations in the (100) direc- 
tion; our choice of units is such that the field, mass and cou- 
pling constants are dimensionless. The term proportional to 
Fa.p(%) = (d — 2)(x 2 8 a — dx a XfAJx d+2 represents the 
dipole interaction. In momentum spacea 

F a(3 (q) = f d d xF a!3 {x)e l i- x = (r 0a + f a q 2 a ) S af3 (3) 
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where ro a , f a , g a p, and h a p depend on details of the under- 
lying lattice. We assume that the non-local quadratic terms 
represented by g a p and h a p obey the same symmetry as the 
local quartic interaction terms u + v a S a p. Thus, in general, 
we will have g a p = g + gi8 a p\ gi > lowers the symmetry 
to Ising. The term ro a combined with local bare mass terms 
makes up r a p in Eq. (|]). We shall consider cubic and tetrag- 
onal symmetry, so r a p — r a 8 a p. In Eq. (^]) and in all of the 
following we use dimensionless momenta q a E [— tt, it]. 

Diagonalization of the quadratic part of the action yields 
the phonon modes, and the paraelectric - ferroelectric transi- 
tion occurs when the lowest zone center mode frequency van- 
ishes. The gradient term in Eq. (Q), along with f a and h, 
controls the dispersion of modes. Note that f a ^ implies 
an anisotropic derivative 5Z Q (V ' a<pa) 2 \ in a spherically sym- 
metric system f a — 0. For simplicity we refer to the case 
gj = 0, f a > as Heisenberg also, because the order pa- 
rameter exhibits a continuous rotational symmetry. Previous 
renormalization group studies of the classical paraelectric - 
ferroelectric-transition have treated the / Q -tejEms as a small 
perturbationEH. Khmel'nitskii and ShneersonEI argue that al- 
though f a in typical materials (e.g. BaTiOa) is of the same 
order of magnitude as g a p and h a p, the anisotropy of observ- 
able quantities is usually weak. Because band theory calcula- 
tions indicate that in many ferroelectric systems f a > 1 are 
quite large we present here a treatment valid for any /. 

The u and v a terms represent local anharmonic interac- 
tions. The materials of main interest here have cubic sym- 
metry in which case v a = v. The quartic interaction in Eq. 
(^) (dropping momentum and energy integrals for simplicity) 
becomes 



(4) 



The term proportional to u is rotationally invariant and in- 
sensitive to the polarization orientation, and the sign of the 
second term determines the polarization orientation in the or- 
dered phase. At the mean field level the action Eq. (^) is 
minimized by polarization of magnitude 
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FIG. 1: Generic ferroelectric propagator Eq. (|ic|) in the static limit 
vj = with cubic symmetry and large anisotropy; q z = 0, r — 0, 
/ = 5.0, A = it. 



When v < Eq. (g) is minimized by a polarization along 
(111) with (f) x = <j>y = cj) z = P/d whereas for v > the 
polarization is along (100) with cj> x = (j> y = 0, <fi z = P. The 
values of the quartic interaction are in each case 



S (4) [0] = 




v < 
v > 
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The condensation energy is of order r 2 /u and as v — * the 
energy barrier separating different symmetry-allowed polar- 
ization directions is a factor of order v ju smaller than the con- 
densation energy. The condition for the stability of a quartic 
interaction is the positive definiteness of Eq. (Q) which (in cu- 
bic symmetry, dimensionality d and at the mean field level) 
translates into 



u + v > 
du + v > 
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If these conditions fail, sixth order terms in <fc have to be in- 
cluded and the transition may be first-order. 



III. PROPAGATOR AND MODES 

Eq. (fy and Eq. (|J) define a model for the phase transition 
in a ferroelectric near a quantum critical point. In the absence 
of nonlinearities, the Heisenberg order parameter correlation 
function G af} = ((j>a4>p)o 
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FIG. 2: One loop diagrams for the renormalization of (a) the quartic 
interaction u and (b) the masses r a . 
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Here w = 2nir(T / 'E n ) is a dimensionless bosonic Matsubara 
frequency, and in our conventions G is dimensionless. 

The nature of the modes defined by the poles of Eq. 
can be best understood by considering the polarization of the 
ferroelectric fluctuation vector <fi. For every q there are d — 1 
transverse and one longitudinal polarizations, all orthogonal 
to each other. The longitudinal mode is always stiff with 
tz7|| = O(g), and h only enters the dispersion of the longitudi- 
nal mode; both g and h are irrelevant to the critical behavior. 
The remaining d — 1 modes are soft and in the case of cubic 
symmetry have the general dispersion 
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where il q is the set of angles defining the direction of q and 
A\ are lengthy expressions derived from Eq. (j8a|). Eq. (||) 
includes all modes with the convention r± = r, r» = r + 
g. For all q such that the polarization of a transverse mode 
points along a crystal axis the respective dispersion softens 
additionally (A\(Q q ) — 0) if / > 0. The effect of / is most 
easily seen by setting q z = and considering only the XY 
block. The resulting transverse mode propagator 



G(m,q) 
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is shown in Fig. ([I]). The dipolar anisotropy (/) leads to ridges 
suggestive of quasi one dimensional behavior. 
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BaTi0 3 


6.36 


5.00 


4.7 


1.25 


0.68 


PbTiOs 


6.79 


5.37 


1.1 


0.26 


-0.09 



TABLE I: Numerical values of the parameters used in the action Eq. 
(^). The values for u and v are initial conditions for the RG flow of 
the two interaction constants in Eq. (Eh. 



IV. INTERACTION RENORMALIZATION 

We obtain the parameters c, r a , f a , g and h in Eq. (||) 
by fitting the poles of Eq. ( |8a| ) to first-principles phonon dis- 
persion curves such as those in RefM We fit the numeri- 
cally calculated mode frequencies near the zone center to the 



FIG. 3: Renormalization diagrams to tree and one-loop order for the 
anisotropic interaction parameter u a p = u + v5 a p in Eq. (j|). This 
is the diagrammatic representation of Eq. (|l l[). 



modes predicted by Eq. (|Ba|) along crystal symmetry direc- 
tions. The soft modes' speed-|C and mass r(T — 0) at the 
lattice constants used in Ref.t3 (ambient pressure) are read- 
ily obtained by fitting the dispersions along (100) to u> = 

\2. 



c^Jr/a 2 + (1 /2)cy/ r / a 2 (q/ a) 2 + 0(q 4 ). The anisotropy pa- 
rameter / is obtained from direct ratios of the curvature of the 
dispersion along (110) and (111). 

The size of the interaction constants u and v can be es- 
timated from first-principles variational studies of a Lan- 
dau free energy of the systemtil E(w) = kw 2 + a'w 4 + 
l' E a >/3 w2 a w % where E(w) is the free energy per unit cell 
and w is a soft mode lattice displacement. The parameters a' 
and 7' (Table V inlill) studied by these authors are related to 
u and v a = v in Eq. (||) by u = (a' + 7' /2)(r / ' n) 2 E and 
v = -(y/2)(r/n) 2 E . ForBaTi0 3 andPbTi0 3 (which have 
cubic symmetry so u a = v and f a — /) we find the results 
listed in Table Q. 

We now study the relevance of quartic interactions in the 
vicinity of the critical point. The one-loop correction to u and 
v is given generically by the diagram in Fig. (||a) and the 
respective renormalization equationEl is given by 

a/3 a/3 

+ 4u; 7(375 0>'V + tWapysftPpt*)] (11) 

where u a p = u + vS a p, the external momentum integrations 
are omitted for brevity, and w Q /3 7< 5 are the one-loop integrals 
over fast modes 
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In Eq. (11)0 = b 1+d / 2 is the field renormalization when the 
fast modes in the shell A/fb<<7<Aare integrated out. The 
diagrammatic version of Eq. ( p"l| ) is shown in Fig. (Jv|). In 
d = 3 the system is in its marginal dimension and the pref- 
actor in Eq. ( |Tl) ) is C / 4 6 _3 ( d + z ) so that the leading interaction 
renormalization is quadratic. The generic form of the renor- 
malization equations is 

dm r d d q 9 , N ru 
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where G is a soft eigenmode of the Gaussian ferroelectric 
propagator. We first consider the simplest case of isotropic in- 
teractions (v a = 0) in an isotropic medium (f a = 0, r a = r) 
in the low temperature limit, and we also let h = for simplic- 
ity. As explained above, the correlation function Eq. (^) then 
has a Heisenberg-like rotationally invariant form with d — 1 
soft eigenmodes G _1 = r + w 2 + q 2 = r + Q 2 and one stiff 
(non-critical) eigenmode G _1 = r + Q 2 + g. Including only 
the soft eigenmodes in Eq. ( pj[ ) the recursion relation for u is 
respectively 
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Eq. (|14j) shows that upper critical dimension of a Heisenberg 
isotropic quantum critical ferroelectric system is d c = 3. 

For a uniaxial (Ising-like) ferroelectric there is a preferred 
'easy axis' for the orientation of <fi which brings about a fur- 
ther increase in effective dimensionality!! The respective cor- 
relation function has the form G{w, q) = (r + w 2 + q 2 + 
QiqI/q 2 ) 1 which gives for Su 
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Thus the upper critical dimension of an Ising isotropic quan- 
tum critical ferroelectric system is reduced to d c = 2. In 
the case of a preferred 'easy plane' of polarization the prop- 
agator has the two eigenmodes G = (r + w 2 + q 2 )^ 1 and 
G = (r + zu 2 + q 2 + g sin 2 6») _1 so that the XY -model and 
Heisenberg ferroelectrics have identical coupling renormal- 
izations which is due to the existence of the same soft mode 
in both cases. 

We now study the possibility that the quasi one dimensional 
behavior associated with / > > 1 may modify the criticality. 
We illustrate the issues using the notationally simpler d = 2, 
g — > oo case, and have verified that our results hold in d = 3 
also. From Fig. (||a) and Eq. ( |lo| ) we obtain, after integration 
over vj and the magnitude of q 



FIG. 4: Coefficients of the angular and radial part of the interaction 
renormalization coefficient Eq. ( |l9[ ) at T = 0. The stable roots of 
B(9) which determine the nature of the fixed points are marked by 
vertical dotted lines. 



the model Eq. (||) should be qualitatively correct except in the 
case of a d — 2 XY ferroelectric, and we exclude this case 
henceforth. 

We further study the fixed points of Eq. ( pi] ) in its full 
anisotropic form. The Gaussian propagator Eq. ( |8a| ) in the 
strong dipole interaction limit g — > oo is 
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In the second, approximate equality we have taken the large / 
limit. We see from this that the quasi one dimensional struc- 
ture does not affect the degree of divergence as r — > 0; in- 
deed / only affects prefactors and not the scale A 2 to which r 
should be compared. To summarize, a mean-field treatment of 
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where Q is defined in Eq. (|b|). With the use of cubic symme- 
try, the possible combinations of iu Q /3 7< 5 are reduced to 

Waf3~/8 = [(Ai - A 2 ) 5 al + A 2 ] S a f3S 7 S 

+ A 3 (1 - S a p) (Sa^Sps + S a sSpy) (18) 

The /- and T-dependent integrals Ax 2,3 are calculated in the 
Appendix. Substituting Eq. (|8j) in Eq. ( [TT1 ) yields coupled 
nonlinear renormalization equations for u and v, similar to 
those written by Aharony and FisherH for the classical case 
(but note that Aharony and Fisher expanded the coefficients 
Ai about the limit of small anisotropy whereas we retain their 
full / dependence). The stability of the Gaussian fixed point 
u = v = is most transparently analysed using polar coor- 
dinates in the (u, v) plane: u — pcosQ and v = psin#. The 
renormalization equations for u and v then become 
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FIG. 5: Phase portrait for Eq. ( |19| ) with / = 1.1. The heavy dashed 
lines indicate the u/v ratios corresponding to stable fixed points of 
Eq. (|l9|); their continuations (light dashed lines) mark the boundaries 
of the unstable region where p flows to large values. The dotted 
lines mark the mean-field stability boundaries Eq. (Q). The shaded 
region indicates initial conditions (uo, «o) which start with Ising-like 
polarization orientation but eventually flow to the fixed line 6\ with 
nearly Heisenberg polarization orientation along (111). 



The 9- and T-dependent coefficients A and B are given in 
the Appendix (Eq. (A7)), and their T — limit is plotted in 
Fig. (Q). These coefficients are to be evaluated at the run- 



ning temperature T(A) = T , p h ys e lnA and are derived on the 
assumption that the physics is dominated by the (Gaussian) 
quantum critical point; in other words, on the assumption that 
control parameter r, interaction amplitude p and temperature 
are not too large. In particular, the model exhibits a pase tran- 
sition at a temperature T c (r) discussed in detail below. At 
temperatures sufficiently near to T c (r) a crossover to physics 
controlled by a classical, non-Gaussian critical point will oc- 
cur and the theory used here ceases to apply. To estimate 
the region of applicability of the equations presented here we 
followB, noting first that the breakdown of the quantum crit- 
ical theory will occur in the classical region T(A) > A. In 
this regime the relevant dimensionless interaction amplitude 
is ^classical = T(A)p(A)/A and Eq. (19a) predicts the classi- 
cal fixed point p* lassical = rim T ^ao Tj[AA{6,T)). We find 
Classical = {1-25318,1.47604} for the two fixed lines 9\, 
6*2 respectively, shown in Fig. (|j), and for / = 1.1. The 
/-dependence of /0*i ass i ca i can be be summarized by the lin- 
ear fits p* lassicl ~ 1.16 + 0.11/ for the 9\ fixed line, and 
Pciassid « 0.48 + 0.77/ for the 9% fixed line. As an estimate of 
the range of validity of the scaling equations, we argue they 
apply for p c i ass icai ^ (l/ 2 )/0* lassical , which corresponds to the 

1/2 

3D Ginzburg criterion p c iassicai/»~ classical ~ E 

We wish to study the stability of the fixed point p = in 
Eq. ([19]). The solutions of Eq. ( |l9| ) asymptotically approach 
the fixed point along 'invariant lines' 9 — given by the sta- 



ble roots 9* of B(6,T) (Eq. (jATbh). Respectively, Eq. 
has a stable fixed point p — > if A(9* , T) > 0. The functions 
A{9) and B{9) are shown in Fig. (|j) for several values of the 
anisotropy /. It is seen that there are three 'invariant lines' 
of which 6>3 (the middle in Fig. fl)) is unstable (A{9t l ) < 0). 
The two stable solutions are 9\ <; corresponding to a nearly 
Heisenberg fixed point \u\ » v; u < 0, and 6% > 7r/2 
corresponding to an Ising-like fixed point v >> \u\, u < 0. 
The dependence of 9\ 2 on / is weak and does not change 
the qualitative behavior. The nature of the 'fixed line' solu- 
tions is most clearly seen in Fig. which shows the phase 
portrait of Eq. (jl9|). The stable fixed lines 9* 2 are shown 
by heavy dashed lines. Above and to the right of the light 
dashed lines the flows are stable (p — > 0); below and to the 
left, unstable (p — > oo). The region of stability found in the 
RG analysis is wider than that found in the mean field approx- 
imation (Eq. (Q) shown in Fig. (^) as light dotted lines). The 
physical content of the two fixed lines 9* 2 is different: 9\ cor- 
responds to a nearly isotropic system with polarization along 
(111) but a relatively weak barrier against polarization reori- 
entation (w/|m| ~ 0.15, but weakly /-dependent), whereas 9\ 
corresponds to a strongly anisotropic system with polarization 
along (100) and a barrier of relative order unity. It is seen in 
Fig. (Q) that there exists a range of initial conditions in the 
shaded wedge between the v = axis and the separatrix in 
the first quadrant, which start with initial values vo > fa- 
voring Ising symmetry but eventually flow to the 9\ fixed line 
with v < and Heisenberg polarization symmetry. 

We see that the ratio and even the sign of u/v may change 
under renormalization. In particular, for initial u > 0, v > 
and u/v less than an (/-dependent) critical value of the order 
of unity, the sign of v changes under renormalization, corre- 
sponding to a predicted change in the polarization direction 
as criticality is approached. Unfortunately, the logarithmic 
nature of the scaling, combined with the n umer ically small 
value of B(9, T) and the factor of p in Eq. (19b) means that 



at T = one must approach criticality extraordinarily closely 
to observe the effect. The scaling turns out to be more rapid 
in the classical regime T(A) > A, but as noted above our 
analysis cannot be extended too far into this regime before the 
equations break down. 

To further illustrate this point and to study how "soon" in 
RG time this change of ordered state orientation occurs, and 
each fixed line is reached, we show in Fig. ^ the evolution 
of the ratio of the two interaction constants along typical tra- 
jectories in Fig. (||). It is seen that the trajectories reach their 
'fixed line' regime relatively late with long temperature de- 
pendent transients sensitive to initial conditions. Trajectories 
that start with Ising symmetry vq > and ultimately flow to 
a Heisenberg fixed line with v < are marked with filled 
symbols. For example, the initial conditions for BaTiC>3 are 
within the shaded range in Fig. (|J), and the sign reversal is 
expected to occur for T = 0.0001i?o at In A — 15. 

The phase portrait Fig. (Q) changes only quantitatively with 
increased anisotropy, e.g. / = 5: the closed loop trajectories 
shrink towards the origin as higher anisotropy reduces both 
interactions u and v (Eq. (\3)), and the slope of the fixed 



lines changes according to the 0-roots shown in Fig. 



6 




CHD ( 0.20,0.20) 
M( 0.01 ,0.02) 
(-0.10,0.10) 
T=0 

T=0.0001 E 
T=0.001E 
T=0.01E„ 



CK3 (0.25,-0.23) 
o—o (0.02,-0.01) 
♦-♦(1.25, 0.68) 
T=0 

— T=0.0001E 
-- T=0.001E 

- T=0.01E„ 



FIG. 6: The solutions of Eq. ( |l9| ) for / = 1.1 plotted for several dif- 
ferent initial conditions. Selected trajectories are shown for several 
temperatures. The trajectories converge into one of two 'fixed lines' 
also seen in Fig. (||). The top and bottom panels show trajectories 
converging to 6*5 and 6*5 respectively. The two panels are shown with 
inverted interaction constant ratios relative to each other to capture 
the significant features of the fixed line approach in each case. The 
arrows mark the quantum - classical crossover ^classical ~ /^classical' 
Trajectories that undergo polarization reorientation are marked with 
filled symbols. 



V. FREE ENERGY, SPECIFIC HEAT AND MASS 
RENORMALIZATION 

Within the Gaussian approximation the free energy per unit 
cell of the system is given by 



(20) 



(21) 



F = Y,I \E w x {q)+T\n(l 

where w\ are the poles of Eq. (p^). The specific heat can be 
obtained directly from this expression as 

c = -t— = v f ddq E^wg^l! 

5T 2 (2n) d smh 2 [zu^q)E /2T} 

Using the general form of the eigenmodes in a cubic Heisen- 
berg system with anisotropy Eq. (m and the isotropic Ising 
mode w 2 (q) = r + q 2 + gi cos 2, 6 the asymptotic low- 
temperature behavior of the specific heat is 

dfi, (T/E ) 3 r°° x 2 (k + x 2 ) dx 
^ [1 + fA x (n q )} 3/2 Jo sinh 2 Vk + x 2 
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FIG. 7: Main figure: critical point specific heat per unit cell to Gaus- 
sian order in units of fcs for a d = 3 Heisenberg model for a set 
of anisotropy parameters /. The vertical arrows mark the points 
t = \J2 + f for each respective curve. Inset: a) the r-dependence of 
C(t) for / = 5; b) Heisenberg to Ising crossover in specific heat, for 
/ = 0. 



T 



t(K + t) 
sinh 2 y/Jt~+t 



dt 



where H and / refer to Heisenberg and Isinj 
and k — (Eo/2T) 2 r. As is seen from Eq. (K 



(22) 

respectively 
) in the low- 
temperature limit the anisotropy / enters only as a multiplica- 
tive factor in a / = expression for the specific heat. The 
specific heats of a d — 3 Heisenberg model are shown in 
Fig. ([7]) for r = (main figure) and away from the critical 
point (inset). We see that except in the unrealistically strong 
(/ > 100) case the only crossover visible is from the quantal 
(C ~ T 3 ) to classical (C ~ const) behavior as T is increased 
through the largest zone boundary phonon frequency (shown 
by arrows in Fig. ([7]). The crossover from Heisenberg to Ising 
symmetry is shown in inset b) of Fig. (J^) for a set of Ising 
interaction strengths gi . 

Finally we study the pressure dependence of the transition 
temperature. The mass flow equation is given by 



dr = 2r(A) dlnA + di?(A) 



(23) 



where di?(A) represents the one loop mass correction Fig. 
(^b). It is possible to express this diagram in terms of an in- 
variant of Eq. (S3) 



dR(A) = -[(d + 2)u(A) + 



(24) 

where the ^-integration and Matsubara frequency summation 
are performed in narrow shells of width dA for each variable 
while the other one is held fixed at the bandwidth cutoff, e.g. 
w n e [A, A + dA] while q = A, and q e [A, A + dA] while 
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FIG. 8: Pressure dependence of the ferroelectric critical temperature 
for the following cases: an Ising model with or = 10 and «o = 1; 
PbTiOs and BaTiOs with parameters in Table U; an isotropic model 
with / = 0, u = 2. 



w n = A. T he tr ace in Eq. (p_4|) is over the propagator eigen- 
values Eq. (Alb). Using the identity 



i 



f°° du> , uj . , , . , 

J ^^th-S^-^q)) 



we obtain the solution of Eq. ([Z3]) for d = 3 in explicit form 



-in \ J / lnA dlnA'e- 21nA '4 



r(A) =e zlnA <>r + 



-u(A')+v(A') 



AVI + Ma^) cothA/2Te lnA ' 



(25) 



coth - 



2Te lnA ' 



l + fA x ((l) 



Here r is the initial condition for the mass; u(A) and u(A) 
are the solutions to Eq. (19); A\(SY) are the angular depen- 
dent anisotropic factors of the dispersion from Eq. (Q), and 
T(A) = Te lnA is the flowing temperature while T is the real 
physical temperature. We find T c from the requirement that at 
the critical temperature the mass flows to zero: r(A — > oo). 
The quantum critical control parameter (which in experimen- 
tal realizations corresponds to e.g. hydrostatic pressure, dop- 
ing, etc.) is r — r c = ra(T) — ro(T = 0). We point out 
that since a 3D ferroelectric is above its upper critical dimen- 
sion we obtain a qualitatively identical phase boundary if we 
simply use the initial conditions u and Vq for the interaction 
constants in Eq. (p5|). In that case T c (r — r c ) can be obtain 
from the simpler expression 



E 

A 



d d q 



[(d + 2)u Q + 3v Q ] 



- 1 



The phase boundary T c (r — r c ) obtained from Eq. (26) is 
shown in Fig. (^) for four representative sets of parameters. 
All curves in Fig. (^J) except Ising behave as T c ~ \r — r^l 2 
near r c ; the Ising behavior is T c ~ \r— r c \ 1/>3 . All curves cross 
over to T c ~ | r — r c \ as T is increased through the softest zone 
boundary phonon frequency Eq . 

We briefly consider the effect of a small density of free car- 
riers characterized by an inter-carrier spacing L p and a dif- 
fusion constant Dp. At length scales longer than Lp and 
frequency scales lower than cup = Dp/L F , these carriers 
will screen the interaction on the scale Lp and overdamp the 
dynamics. The details of the crossover depend on the ra- 
tio Dp/cLp. If Dp/cLp >> 1, then there is a two-stage 
crossover: as the scale is decreased, first the dynamics be- 
comes overdamped and then subsequently the characteristic 
length scale passes through the screening length and the inter- 
action becomes effectively short ranged. On the other hand, if 
Dp/cLp << 1, then screening and o verdamping occur at the 
same scale. Further studies of this crossover will be presented 
elsewhere. 

All cases except for 2d XY symmetry are above the upper 
critical dimension enabling a controlled treatment. Lattice- 
induced anisotropies arising from the dipolar interaction are 
not small in real materials, and lead e.g. to strong "quasi 
one-dimensional" effects in the phonon spectrum (cf Fig. [I]). 
However, we showed that for systems above the upper criti- 
cal dimension the effect on the critical behavior is unimpor- 
tant; only for unrealistically strong anisotropies / > 100 is 
an intermediate quasi one-dimensional regime visible in the 
specific heat. A change of polarization direction under scal- 
ing is suggested for BaTiOs near T c (p) and for p sufficiently 
close to the critical pressure at which T c — > (although the 
scaling equations break down at approximately the scale of 
the anisotropy change). We have presented exact results, in 
physical units, for the phase boundary and specific heat. For 
PbTi03 and BaTiO,3 quantum critical effects are dominant for 
T < 50K if the materials are tuned by pressure to the quantum 
critical point. 
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VI. APPENDIX 

We calculate the one-loop diagrams Eq. (|l2|) by using a 
diagonal representation for the Gaussian propagator achieved 
through the rotation matrix R 



Gap - Raag<j{^n,q)R a l 

where 

g a {m n ,q) 



(Ala) 



1 



1 



wl+u)l{q) wl+r a + q 2 (l + fA a (9,ip)) 

(Alb) 



(26) The one-loop integrals then become 
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i 



-R: 



/-dw /■ d d g v x x 

u x 7 /if L 



i 



(A2) 



coth • 



w a +w*(g)][w a +a;2(g)] 2T 



We perform the integration over the magnitude of q in d — 3 and obtain the remaining integrals over angles only, which we then 
calculate numerically 



16^3 A E J ' d^R^R^^Rj f[Av( pj_ A 



-AM] 



A 



2T 



coth — coth ; 



AVTT7A; { cQth A \ f 0r am ± A v (n q ) 



2T 



2T 



1 dA 
16tt 3 X 



2[i + M M (n,)r 



C0 thM±Z5 M 

2T 



sinh^/1 + /.1,,(<2„) 



, A 
coth — 
2T 



The isotropic case / = is described by the second expression in Eq. (43). For numerical calculations Eq. ([\3|) is more 
conveniently written as 



for A^Slq) = A v {VL q ) 



(A3) 



A f 



Wa/3-rS = < 



.9m + 

^cothA/(2T^) - ^ coth A/ (22^) 



Y^jdA^ y dflqRa^Rp^R^vRsv ^~ 



, A 
coth — 
2T 



4 /2 



coth 



A 



1 



A/(r^) 

2T^ ^ ' sinhA/(T^) 
I 



coth 



A_ 

2T 



for 5 P = 



(A4) 



Here R afi is a matrix whose columns are the fi-th eigenvector 
and — A 2 [l + fA^(fl q )] is the /^-th eigenvalue of Eq. ( |l7| ) 
both evaluated at tu = 0; r = 0; q = A, and both having an 
i mpl icit angular dependence. At low temperatures T < A Eq. 
( A4 ) reduces to 



^-jdA^ [ &VL q R atl Rp ll R lu R SlJ 

\1V 



V9^ 
(A5) 

The cubic-symmetric integrals Ai.2,3 appearing in Eq. ( |T8| ) 
are defined as 



A 2 (b)=T^J-^ d G™Gyy = w xxyy 
Mb) =Tj2j-0^ (G xy ) 2 = w xyxy (A6c) 



(A6a) 
(A6b) 



and are calculated from Eq. ( A4) numerically. 

The angle-dependence of the p and 9 derivatives in Eq. ( |l9| ) 
is given by 

A(9, T) = ax cos 3 6 + c 2 sin 3 9 (Ala) 
+ sin 9 cos [(61 + 02) cos 6* + (ci + 62) sin#] 

B(0, T) = a 2 cos 3 - c x sin 3 6» (A7b) 
+ sin#cos#[(&2 — cos9 + (c 2 — 61) sin#] 

where 

01 = 4(3,4! + 4A 2 + 144a) a 2 = 16(Ai - A 2 - 2A 3 ) 
61 = 8(5Ai + 16A 3 ) 6 2 - 48(Ai - A 3 ) 

ci = 36A 3 c 2 - 36(Ai - A3) 



9 



Lines M.E. and Glass A.M., Principles and Applications ofFerro- 
electrics and Related Materials (Oxford University Press, Oxford 
1977) 

K. A. Muller and H. Burkard, Phys. Rev. B 19, 3593 (1979) 
Hertz, JA. Phys. Rev. B 14, 1165 (1973) 

Sachdev, S. Quantum Phase Transitions (Cambridge University 
Press 1999)) 

Millis, A.J. Phys. Rev. B 48, 7183 (1993) 

A. B. Rechester, Sov. Phys. JETP 33, 423 (1971) 

D. E. Khmel' nitskii and V. L. Shneerson, Sov. Phys. - Solid State 

13, 687 (1971) 

T. Moriya, Spin Fluctuations in Itinerant Electron Magnets 



(Springer- Verlag, Berlin 1985) 

9 Amnon Aharony and Michael E. Fisher, Phys. Rev. B 8, 3323 
(1973) 

10 Ph. Ghosez, E. Cockayne, U. V. Waghmare, and K. M. Rabe, 
Phys. Rev. B 60, 836 (1999) 

11 R. D. King-Smith and David Vanderbilt, Phys. Rev. B 49, 5828 
(1994) 

12 A. Aharony, Phys. Rev. B 8, 3358 (1973) 

13 R. Yu and H. Krakauer, Phys. Rev. Lett. 74, 4067 (1995) 

14 Waghmare U.V. and Rabe K.M., Phys. Rev. B 55, 6161 (1997) 



